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Abstract. Workspace and joint space analysis are essential steps in describing the task and de¬ 
signing the control loop of the robot, respectively. This paper presents the descriptive analysis of 
a family of delta-like parallel robots by using algebraic tools to induce an estimation about the 
complexity in representing the singularities in the workspace and the joint space. A Grobner based 
elimination is used to compute the singularities of the manipulator and a Cylindrical Algebraic 
Decomposition algorithm is used to study the workspace and the joint space. From these alge¬ 
braic objects, we propose some certified three dimensional plotting describing the the shape of 
workspace and of the joint space which will help the engineers or researchers to decide the most 
suited configuration of the manipulator they should use for a given task. Also, the different pa¬ 
rameters associated with the complexity of the serial and parallel singularities are tabulated, which 
further enhance the selection of the different configuration of the manipulator by comparing the 
complexity of the singularity equations. 
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1 Introduction 

The workspace can be defined as the volume of space or the complete set of poses 
which the end-effector of the manipulator can reach. Many researchers published 
several works on the problem of computing these complete sets for robot kinemat¬ 
ics. Based on the early studies |[2T][T2l, several methods for workspace determi¬ 
nation have been proposed, but many of them are applicable only for a particular 
class of robots. The workspace of parallel robots mainly depends on the actuated 
joint variables, the range of motion of the joints and the mechanical interferences 
between the bodies of mechanism. There are different techniques based on geomet¬ 
ric HU El, discretization llll|3l[9l, and algebraic methods [17] [5l [3 which can 
be used to compute the workspace of parallel robot. The main advantage of the ge¬ 
ometric approach is that, it establish the nature of the boundary of the workspace 
113 . Also it allows to compute the surface and volume of the workspace while being 
very efficient in terms of storage space, but when the rotational motion is included. 
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it becomes less efficient. 

Interval analysis based methods can be used to compute the workspace but the com¬ 
putation time depends on the complexity of the robot and the requested accuracy 
191 . Discretization methods are usually less complicated and can easily take into ac¬ 
count all kinematic constraints, but they require more space and computation time 
for higher resolutions. The majority of numerical methods which is used to deter¬ 
mine the workspace of parallel manipulators includes the discretization of the pose 
parameters for computing workspace boundaries lH. There are other approaches, 
which are based on optimization algorithms 1 ^ for fully serial or parallel manipu¬ 
lators, analytic methods for symmetrical spherical mechanisms O. In (Tl a method 
for computing the workspace boundary for manipulators with a general structure 
is proposed, which uses a branch-and-prune technique to isolate a set of output 
singularities, and then classifies the points on such set according to whether they 
correspond to motion impediments in the workspace. A cylindrical algebraic de¬ 
composition (CAD) based method is illustrated in Eia, which is used to model the 
workspace and joint space for the 3 RPS parallel robot. 

This paper presents the results which are obtained by applying algebraic methods 
for the workspace and joint space analysis of a family of delta-like robot including 
complexity information for representing the singularities in the workspace and the 
joint space. The CAD algorithm is used to study the workspace and joint space, 
and a Grobner based elimination process is used to compute the parallel and serial 
singularities of the manipulator. The structure of the paper is as follows. Section 2 
describes the architecture of the manipulator, including kinematic equation and joint 
constraints associated with the manipulators. Section 3 discusses the computation of 
parallel as well as serial singularities and their projections in workspace and joint- 
space. Section 4 presents a comparative study on the shape of the workspace of 
different delta-like robots. Section 5 finally concludes the paper. 


2 Manipulators Architecture 


The manipulator architecture under the study is a three degree of freedom paral¬ 
lel mechanism which consists of three identical legs, the different arrangements of 
these legs give rise to family of delta like robot. Several types of delta-like robot 
were studied, few of them are Orthoglide Elm, Hybridglide, Triaglide ca and 
UraneSX (91. The position vectors of end points of leg are and B^, also Aj 
and are the position vectors of end points of actuator where j represents the 
manipulator type from the family of delta like robot (7 = 1 - Orthoglide, 2 - Hy¬ 
bridglide, 3 - Triaglide, 4 - UraneSX). pj represents the prismatic joint variables 
whereas represents the position vector of the tool center point which is shown in 
Eq. 0. 

||A/b/||=p/ p=[xyzf with /= 1,2,3 and 7 = 1,2,3,4 (1) 
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Fig. 1 Configuration plot for Orthoglide (a), Hybridglide (b), Triaglide (c) and UraneSX (d) robot 


The kinematic equations of the family of delta like robot can be generalized as 
||P —B/II = L/. All the computations and analysis are done for L/ = L and by 
imposing the following constraints on joint variables. Without joint limits, all the 
robots admit two assembly modes and eight working modes. 

0 < Pi < 2L 0 < p2 < 2L 0 < P3 < 2L (2) 

The Ortho glide mechanism is driven by three actuated orthogonal prismatic joints. 
A simpler virtual model can be defined for the Orthoglide, which consists of three 
bar links connected by the revolute joints to the tool center point on one side and 
to the corresponding prismatic joint at another side. Several assembly modes of 
these robots depends upon the solutions of direct kinematic problem is shown in 
Fig. [Ja). The point represents the pose of corresponding robot. However more 
than one value of i for the point P/ shows the multiple solutions for the DKP. The 
constraint equations for the Orthoglide are 

{x- Pi)'^+y^+z^ =L^ ■,x^+{y-p2f+z^=L^-,x^+/+{z-p3f=L^ 

The Hybridglide mechanism consists of three actuated prismatic joints, in which 
two actuators are placed parallel and third one perpendicular to others two. Also 
the three bar links connected by spherical joints to the tool center point on one side 
and to the corresponding prismatic joint at another side. Several assembly modes 
of these robots depends upon the solutions of direct kinematic problem is shown in 
Fig.[2b) and following is the constraint equations: 
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(x-lf + iy-pif + z^=L\ix+lf + iy- P2f+z^ = lW +/ + (^ - Ps)" = 

The Triaglide manipulator is driven by three actuated prismatic joints, in which all 
the three actuators parallel to each other and placed in the same plane. The architec¬ 
ture of the Triaglide is shown in Fig. [^c) and the constraint equations are defined 
as: 

{x-\Y + {y- Pif + £-=L^\{x+\f + {y- p2f + + {y- p^f+ 7 ^=L^ 

The UraneSX is similar to triaglide, but instead of three actuators in the same plane, 
it is placed in different planes. The architecture of the Triaglide is shown in Fig.[2d) 
and the constraint equations are as follows: 

(v — 1)^ + (z — Pi)^ =l}\{x-\r 1/2)^ + (y ~ + (z — =1^\ 

(x+ 1/2)^ + (y + a/3/2)^ + {z — Ps)^ =L^ 


3 Singularity Analysis 


Singularities of a robotic manipulator are important feature that essentially influence 
its capabilities. Mathematically, a singular configuration may be defined a rank de¬ 
ficiency of the Jacobian describing the differential mapping from the joint space to 
the workspace and vice versa. Differentiating the constraint equations of the robot 
with respect to time leads to the velocity model: At + Bq = 0 where A and B are the 
parallel and serial Jacobian matrices, respectively, t is the velocity of P and q de¬ 
fines the joint velocities. The parallel singularities occur whenever det(A) = 0 and 
the serial singularities occur whenever det(B) = 0 |[ 8 l. 

det(A)o = -8pip2p3 + 8pip2Z + 8pip3y + 8p2p3X 

det(A)/, = - 8 pip 3 V + 8 p 2 p 3 X- 8 pip 3 + 8 piz - 8 p 2 p 3 + 8 p 2 Z+ 16p3y 

det(A)^ = 8 piz + 8 p 2 Z- 16p3Z 

det(A)^ = 4\/3(3z - pi - p2 - p3 + psx^pix - 2pix) -f 12p3y - 12p2y (3) 

Parallel and serial singularities as well as their projections in workspace and joint 
space are computed using a Grobner based elimination method. This usual way for 
eliminating variables (see fTOl ) computes (the algebraic closure of) the projection 
of the parallel singularities in the workspace. In the same way, one can compute 
(the algebraic closure of) the projection of the parallel singularities in the joint 
space. Both are then defined as the zero set of some system of algebraic equa¬ 
tions and we assume that the considered robots are generic enough so that both 
are hypersurfaces. det(A)c,, det(A)/j, det(A)^ and det(A)j^ are the parallel singulari¬ 
ties of Orthoglide, Hybridglide, Triaglide and UraneSX, respectively. Starting from 
the constraint equations and the determinant of the Jacobian matrix, we are able to 
eliminate the joint values. This elimination strategy is more efficient than a cascad- 
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ing elimination by means of resultants which might introduce many more spurious 
solutions : singular points that are not projections of singular points. Due to the lack 
of space, other equations associated with serial singularities are not presented in this 
paper. Figurej^represents the projections of parallel singularities in projection space 



Fig. 2 Projection of parallel singularities of Orthoglide (a), Hybridglide (b), Triaglide (c) and 
UraneSX (d) in workspace. All the computations are done with joint limits. 


(x,y,z). These surfaces represent the singularity associated with the eight working 
modes. 


Table 1 Comparison of the different parameters associated with the singularities for the robots 


Manipulators Types of 

Singularities 

Plotting Time 
(s) 

Degrees 

No.of terms Binary 
Size 

No. of Cells 

Orthoglide 

Parallel 

0037.827 

18[10,10,10] 

097 

015 

[02382,0272] 


Serial 

0005.133 

18[12,12,12] 

062 

012 

[00044,0004] 

Hybridglide 

Parallel 

5698.601 

20[16,08,12] 

119 

017 

[28012,1208] 


Serial 

0007.007 

18[12,12,12] 

281 

017 

[00158,0027] 

Triaglide 

Parallel 

0010.625 

03 [00,00,03] 

002 

002 

[00138,0004] 


Serial 

0005.079 

06[06,06,06] 

042 

007 

[00077,0017] 

UraneSX 

Parallel 

0022.625 

06[06,04,00] 

015 

040 

[02795,0070] 


Serial 

0018.391 

12[12,12,12] 

252 

151 

[00392,0142] 


In Table a comparative study of five parameters among the family of delta like 
robot is presented. We have tabulated the main characteristics of the polynomials (In 
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three variables) used for the plots (Implicit surface): their total degree, their number 
of terms and the maximum bitsize of their coefficients. We have also reported the 
time (In seconds) for plotting the implicit surface which they define and the number 
of cells computed by the CAD, as well as the number of cells in the final result after 
gluing those that are adjacent and belongs to the same connected component. Sev¬ 
eral functions are used which involves the discriminant variety, Grobner bases and 
CAD computations, computed in Maple 18 with a Intel(R) Core(TM) i7-3770 CPU 
@ 3.40 GHz (14 Gb RAM). As can be seen from Tablethere exists higher values 
of all the parameters for the Hybridglide, among all manipulators listed, which in¬ 
fers that it has more complex singularities, whereas for the Triaglide all the values 
are least which intuits the less complicated singularities. For example, the compu¬ 
tation times for the Hybridglyde for parallel singularities is high compared to the 
one for the Othoglide, even if the surface has similar characteristics. This is due 
to the geometry of the surface which is more difficult to decompose in the case of 
the Hybridglide : the Cylindrical Algebraic Decomposition is described by 1208 
cylindrical cells in the case of Hybridglide while it is described by 272 cells for the 
Orthoglide. 


4 Workspace analysis 


The workspace analysis allows to characterize of the workspace regions where the 
number of real solutions for the inverse kinematics is constant. A CAD algorithm 
is used to compute the workspace of the robot in the projection space (v,y,z) with 
some joint constraints taken in account. The three main steps involved in the analysis 
are 

• Computation of a subset of the joint space (resp. workspace) where the number 
of solutions changes: the Discriminant Variety . 

• Description of the complementary of the discriminant variety in connected cells: 
the Generic Cylindrical Algebraic Decomposition (CAD). 

• Connecting the cells belonging to the same connected component in the counter¬ 
part of the discriminant variety: interval comparisons. 

The different shapes of workspace for the delta-like robots is shown in Fig.[^ where 
blue, red, yellow and green regions correspond to the one, two, four and eight num¬ 
ber of solutions for the IKR A comparative study is done on the workspace of 
the family of delta-like manipulator and the results are shown in Fig. All the 
workspace are plotted in the rectangular box, where x G [—2,2], y G [—2,6] and 
z G [—2,6], so that the shapes of these workspace can be compared. From the Fig.[^ 
it can be intuited that the Triaglide will be good selection, if the task space is more 
in horizontal plane, whereas the Orthoglide is good for the three dimensional task 
space. 
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(b) 


(d) 


Fig. 3 Workspace with joint constraints for Orthoglide (a), Hybridglide (b), Triaglide (c) and 
UraneSX (d) robot with the number of inverse kinematic solutions. Blue, red, yellow and green 
regions corresponds to the one, two, four and eight number of solutions for the IKR 


5 Conclusions 

A comparative study on the workspace of different delta-like robots gives the idea 
about shape of the workspace, which further plays an important role in the selec¬ 
tion of the manipulator for the specific task or for the trajectory planning. The main 
characteristics associated with the singularities are tabulated in Table which also 
gives some information about the complexity of the singularities, which is an essen¬ 
tial factor for the singularity-free path plannings. From these data, it can be observed 
that the singularities associated with the Hybridglide are complicated, whereas the 
structure of those associated with the Triaglide is rather simple. 
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